Time resolved pattern evolution in a large aperture laser 



F. Encinas-Sanz, I. Leyva and J. M. Guerra 
Dept. Optica, Facultad de CC. Fisicas 
Universidad Complutense de Madrid 28040 Madrid, Spain 



Abstract 



We have measured quasi-instantaneous transverse patterns in a broad 
aperture laser. Non-ordered patterns yielding to boundary determined regular 
structures in progressive time-integrated recording are observed. The linear 
analysis and numerical integration of the full Maxwell-Bloch equations allow 
us to interpret the features of the experiment. We show that this system being 
far from threshold cannot be fully understood with a perturbative model. 
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Pattern formation in systems which exhibit spatio-temporal chaos has been a field of 
intense research in last years. An averaging process from chaotic to boundary-selected or- 
dered patterns has been observed in hydrodynamics . A similar behaviour is foreseeable 
in other pattern forming systems, such as large aperture lasers. 

In fact, this phenomena has been predicted for lasers from the basis of the Maxwell-Bloch 
equations 0-^. However, due to its extremely fast evolution, the time resolved spatio- 
temporal dynamics of a broad area laser has never been observed. So far the experimental 
work only analyze averaged patterns, mainly in cw CO2 0)@H11 semiconductor lasers , 
since the minimum reachable exposure time was about l/xs, too long to obtain information 
about the pattern history. 

In this letter we study, both experimentally and theoretically, the time resolved dynamics 
of a large Fresnel number pulsed laser. The measurements were carried out with a system 
that we have recently developed to obtain infrared snapshots with a short exposure time 
(minimum ~ 1 ns). This setup has been described in detail in Ref. and allows us to 
record virtually instantaneous transverse intensity laser patterns. 

The source was a transversely excited atmospheric (TEA) CO2 laser, with Fresnel number 
jr = |_ ^ 10, where 2b = 20 mm is the laser aperture, L = 1 m the resonator length and 
A = 10.6 X 10~^m the lasing wavelength. This laser emits about 15% of the output energy 
in a gain switch pulse (~ 70 ns), followed by a long collisional transfer tail (~ 2-3 /is) 



TT[ . Being a pulsed laser, it offers interesting characteristics such as high pumping far from 
threshold and wider detuning, which allow the exploration of additional features of the class 
B laser behaviour. 

In order to study the dynamics of the system, we take snapshots all along the duration 
of the laser pulse. A sample of recorded patterns taken with a 6 ns temporal resolution at 
different times along the pulse is shown in Fig. Note that all the instantaneous patterns 
are disordered and are non-reproducible from shot to shot (Fig. |l]a, b, c, (color)). 

However, such irregular appearance masks some kind of regularity, since the averaged 
patterns integrated all along the pulse length are ordered and reproducible (Fig. |l|d), having 
eight or nine rolls parallel to the almost flat laser electrodes. The transverse spatial period 
of those bands is ~ 1.8 mm, similar to the size of the intensity maxima appearing in the 
disordered patterns. No sign of two spatial structure scales can be observed, contrary to 
the predictions of Ref. 0. It is remarkable that the same regular structure is recovered by 
averaging over many equivalent instantaneous patterns, observation that was suggested by 
the results in Ref. |]r2|. In this sense, a sort of ergodicity is observed. 

We have also studied the temporal evolution of the intensity in one small area (~ 1 
mm^), for which a photon-drag detector was used (rise time ~ 1 ns). We flnd that the 
local intensity oscillates in a completely irregular form (Fig. 0a), with a period of 10 ns 
approximately. Furthermore, the cross correlation between the local oscillation measured at 



two different points of the patterns is very low, even if they are as close as 6 mm |T^ . 

As the characteristic period of the intensity local fluctuations is about 10 ns, in order 
to record true instantaneous snapshots it would be convenient to reduce the width of the 
temporal window as much as possible below this value. But, if the window width is too 
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small we do not get enough photons to be recorded by our system. On the other hand, a 
snapshot recorded with a much longer exposure integrate the pattern over several periods 
and it cannot be considered instantaneous any longer. That justify the 6 ns choice that we 
have made for the temporal window width. 

Nevertheless, it would be interesting to know how much recording integration time is 
necessary to obtain a sufficiently ordered pattern. Therefore, in order to follow the averaging 
process, we have also measured sequences of progressively integrated patterns by varying the 
temporal window (Fig. Q). These measurements were made in the first 150 ns of the laser 
pulse, in which the intensity is several times larger than in the pulse tail, and therefore the 
limitation due to the small intensity disappears. That allows us to reduce the time window 
down to 2 ns (Fig. ^), showing a small number of intensity maxima which contrast deeply 
with the background. We see that the longer exposure, the bigger number of maxima appear 
in the (Fig. |]b, c), approaching progressively the ordered pattern, whose rolls are already 
clearly recognizable (Fig. ^d) when the time integration is about 100 ns, yet much shorter 
than the total pulse duration. 

Our theoretical approach to the problem includes both numerical simulations and a linear 
stability analysis of the full Maxwell-Bloch laser equations. 

Hence, in order to reproduce the spatio-temporal dynamics, we directly integrate the 
two- level Maxwell-Bloch equations 0: 

k[{1 - z5) - i^At]E - kQP , (1) 
^4DE+{l + t5)P] , (2) 
lll[D-l-^iEP* + E*P)], (3) 

where E=E(x,t) is the slowly varying electric field, P=P(x,t) the polarization, D=D(x,t) 
the population inversion, Q=Q(x,t) the rescaled pump, k, = ^ln{R) representing the cavity 
losses with R = ^/R^R^ = 0.78 the resonator refiectivity , a = -^^^^ a diffraction coefficient, 
^ _ ^21-^ ^]^g rescaled detuning and the Laplacian in the adimensional transverse coor- 
dinates of the system x=(x,y). In an atmospheric laser, the decay rates can be chosen as 
7^ = 3 X 10^ s-\ 7|| = lO^s-i and /c = 3.9 x 10^ s-^. 

The transversal pumping profile is taken to be homogeneously distributed along one 
transverse axis and Gaussian in the other, in order to reproduce the experimental current 
spatial distribution of the discharge. Likewise, the temporal form of the pumping was 
simulated by a function approximating the pulsed excitation, with typically large pumping 
parameters in the maximum (Qmax —15-25). Null boundary conditions, equivalent to the 
experimental ones, were used. 

The patterns and local intensity temporal evolution obtained by a standard numerical 
method, look very similar to those recorded experimentally. In the spatial domain, the 
instantaneous patterns (Fig. |l|e) are well reproduced, as well as the time integrated (Fig. 
|l]f). In both them, only one spatial scale of structures is found, in accordance with the 
experimental records. The numerical temporal evolution (Fig. Qd), although completely 
irregular, presents a characteristic time of the same order than its experimental counterpart 
(Fig. 0a). 
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As it has been mentioned, this time scale of the irregular oscillations measured in the 
local dynamics is around a few nanoseconds (Fig. whereas the typical evolution time of 
the pulsed pumping is a few microseconds (2-3 /xs). A good approximation is to consider 
the dynamics measured in the slow varying pulse tail to be quasi-stationary |jT^. Thus, we 



can still use the properties of a linear stability analysis to gain some deeper insight into 
the observed irregular dynamics. Taking this consideration into account, we undertake the 
stability study linearizing Eqs.(l-3). As it is known, the Liapunov exponents of the equations 
are the real parts of the five roots of the corresponding secular quintic equation, which cannot 
be solved analitically, and therefore it is not possible to obtain algebraic expressions for the 
eigenvalues. However, it is known that, for class B lasers, out of the five roots, one is real 
and the rest come in complex conjugated pairs 

It can be shown that the resonator can develop the instability associated to the real root 
only when the condition 

^ Q An^) ' 

is satisfied Q. In a TEA-CO2 laser, only one molecular transition oscillates (P20 line, Ref. 



rSp. In the present case, it includes a large number of axial modes simultaneously (15-20), 
with a free spectral range of ^ ~150 MHz [jl3|. Hence, most of the values of S lie in the 



interval —1<S<1. Since Q ^ 1, the condition (^ is not satisfied for the most of the 
oscillating axial modes. Thus, we conclude that this instability does not actually affect the 
laser, or does so very weakly. 

In addition, two of the complex conjugate roots can be approximated as — 7±(1 ± i5), 
being their real part always negative and therefore not associated to any instability. 

The detuning value determines whether the real parts of the two remaining conjugate 
roots are negative for every transverse wave-vector k, or positive for an interval around a 
value ko{S), the wavelength vector at which the Liapunov exponent is maximum and positive 
(Fig. I). Each set of system parameters has a certain critical detuning Sc, such as for 6 > 6c, 
there is an interval of wave numbers k for which the real part of the root is positive. In the 
present case this critical value is rather low {dc ~ 0.06). Then, for the most the axial modes 
with 6 > 0, the instability associated to those complex conjugate roots shows up. 

Summarizing, in this kind of laser most of the axial modes with positive detuning bear 
the short- wavelength instability due to the pair of complex conjugate roots, whereas that 
associated to the real root is not supported because of the diffraction (Fig.^). Therefore, 
the observed irregular spatio-temporal behaviour can only be justified by the action of the 
remaining instability, in contrast with ^j. 

On this basis, it is possible to estimate the expected spatial and temporal scales. The 
spatial scale of the instabilities associated to the complex root will be around k~^. By solving 
numerically the secular equation with d = 0.6 and Q = 6.0 (a reasonable mean value in the 
pulse tail), we obtain kg = 3225 m~^ Then the size of the generated structures should be 

So = 2TxkQ^ ~ 1.94 X 10"^ m (5) 

The average size of the experimentally measured intensity maxima in the instantaneous 
patterns, and consequently the spatial period of the bands appearing in the time integrated 
patterns is 
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Sexp = 1.8 X 10-'^ m (6) 

The agreement between (|]) and (^) relates the observed dynamics with the remaining 
instabihty. 

Furthermore, the imaginary part of these eigenvalues gives the oscillation frequency of 
the solutions. Hence, for the same ko, we obtain a; ~ 360 x 10^s~^, corresponding to a period 
of 16 ns, very close to both the experimental and the numerically obtained fluctuation period 
of ~ 10 ns (Fig. I^a, b). Thus, the time scale is also well predicted by the stability analysis. 

Concerning the physical interpretation of the irregular observed dynamics, in the present 
case the suppression of the instability coming from the real root invalidates its being con- 
sidered as the origin of the disordered instantaneous pattern. Thus we found that even 
though the distance to the threshold is moderate here {Q ~ 6), an analysis based on the 
order parameter equations seem to be already insufficient. To test this, the same numerical 
integration has been carried out for an hypothetical case nearer threshold {Qmax = 2), where 
the order parameter equations must dominate the amplitude behaviour [Q. In agreement 
with this perturvative approach, since the phase instability is inhibited here, the instanta- 
neous patterns show a large degree of order (regular rolls which oscillate periodically. Fig. 

However, in the far threshold pumping case they are disordered, and therefore, the 
present problem does not seem suitable for reduction to a perturbative one, as is usual in 
the theoretical approaches. In other words, in the observed dynamic it is not possible to 
distinguish between phase and amplitude fluctuations. This result is a test of the validity 
range of the order parameter equations, which was not easily verified experimentally. As 
a more probable origin of the phenomenon, a secondary instability of the traveling wave 
solutions can be suggested |T6[ . 



In conclusion, in this work we report the measurement of time resolved intensity patterns 
in a large aperture laser, by means of an experimental system developed in our laboratory. 
A rich irregular intensity spatio-temporal dynamics, usually masked under time integrated 
measurements, has been uncovered. We show experimental evidence of how this local ir- 
regular dynamics averages to boundary determined order, as had been observed in other 
pattern forming systems but so far only predicted theoretically for lasers. Besides, a numer- 
ical integration of the two-level full Maxwell-Bloch equations and its corresponding stability 
analysis reproduce the experimental observations with outstanding agreement. 



5 



FIGURES 



FIG. 1. Instantaneous patterns at different times of the laser pulse. Recordered: a) 150 

ns b)300 ns c)500 ns delay from the gain-switch pulse. Numerically generated: e)500 ns delay 
with Qmax=^8.0 and 6 = 0.6]. Time-integrated patterns [experimental, d), and numerical f)]. 
Experimental pattern dimension 20x20 mm. 

FIG. 2. Time evolution of the local intensity, a) measured b) numerically generated for 
Qmax = l8.0,S = 0.6. 

FIG. 3. Patterns recordered with different exposure time: a)2 ns , b) 6 ns c)30 ns d) 100 ns 

FIG. 4. Real(A) for several 6 values, being the diffraction limit for the instabilities, and 
Kexp the experimentally found mean wave number. 

FIG. 5. Numerically generated instantaneous pattern for Qmax='^-0, 6=0.6. 
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This figure "FIGl.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/patt-sol/9911007vl 



This figure "FIG2.GIF" is available in "GIF" format from: 
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